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Abstract 

Context. The recent findings of gas giant planets around young A-type stars suggest that disks surrounding Herbig Ae/Be stars will 
develop planetary systems. An interesting case is HD 142527, for which previous observations revealed a complex circumstellar 
environment and an unusually high ratio of infrared to stellar luminosity. Its properties differ considerably from other Herbig Ae/Be 
stars. This suggests that the disk surrounding HD 142527 is in an uncommon evolutionary stage. 

Aims. We aim for a better understanding of the geometry and evolutionary status of the circumstellar material around the Herbig 
Ae/Be star HD 142527. 

Methods. We map the composition and spatial distribution of the dust around HD 142527. We analyze SEST and ATCA millimeter 
data, VISIR N and Q-band imaging and spectroscopy. We gather additional relevant data from the literature. We use the radiative 
transfer code MCMax to construct a model of the geometry and density structure of the circumstellar matter, which fits all of the 
observables satisfactorily. 

Results. We find that the disk of HD 142527 has three geometrically distinct components separated by a disk gap running from 30 to 
130 AU. There is a geometrically flat inner disk running from 0.3 AU up to 30 AU; an optically thin halo-like component of dust in 
the inner disk regions; and a massive self-shadowed outer disk running from 130 AU up to 200 AU. We derived a total dust mass in 
small grains of 1.0-10"'' M© and a vertical height of the inner wall of the outer disk of h = 60 AU. Owing to the gray extinction of the 
'halo' we obtained new stellar parameters, including a stellar luminosity of 20+2 Lq and age of lO*'*"* yr 

Conclusions. We find that the disk surrounding HD 142527 is highly evolved despite the relatively young age of the star. The peculiar 
disk geometry can explain the extreme IR reprocessing efficiency of the disk. Furthermore, the geometry, the large disk mass, and the 
highly processed dust composition are indicative of on-going planet formation. 

Key words. Stars: formation - Protoplanetary disks - Stars: variables: T Tauri, Herbig Ae/Be - Radiative transfer - HD 142527 



1. Introduction 

i Herbig Ae/Be (HAeBe) stars are intermediate mass, pre-main 
sequence (PMS) stars with an infrared (IR) excess and emis- 
sion Unes. The observed IR excess is caused by circumstellar 
, dust confined to a disk. Based on the amount of far-IR excess 
' emission the geometry of these disks is considered to be either 
flaring or flattened, classified by Meeus et al. (2001) as group I 
or II respectively. The composition, structure and evolution of 
these disks has been extensively studied because they are gener- 
ally believed to be the sites of ongoing planet-formation. Recent 
discoveries of gas giant planets in wi de orbits around interme- 
diate mass young st ars (iMarois et al.. 200$; Kalas et al, 2QQ8^ 
iLagrange et al.ll2009l) provide exciting support for this interpre- 
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tation of disks around Herbig Ae/Be stars. In this study, we focus 
on the F6IIIe star HD 142527. This star is particularly interesting 
because its extremely large IR and millimeter e xcess is hard to 
understand in terms of a passive disk model (see lDominik et all 
i2003ft . 

The disk of H D 142527 has been imaged in scattered light 
dFukagawa et alJ 12006) and in the thermal IR dFuiiwara et al.l 
2006), showing that the disk extends to several hundred AU with 
a prominent disk gap at a distance of about 100 AU. The dust in 
the disk is highly processed. Both the scattered light and thermal 
IR images show evidence for grain growth: the inferred typical 
grain sizes are 1 -2 fim. This picture is also seen in the thermal 
IR spectrum of HD 142527. Van Boekel et al. (2004a) show that 
the silicates in the inner disk, spatially resolved with the MIDI 
instrument at the Very Large Telescope Interferometer, are fully 
crystalline. Even the integrated disk spectrum of HD 142527 has 
a crystallinity of more than 20%. The outer cold dust is char- 
acterized by crystal line water ice and possibly hydro-silicates 
dMalfait et al.lll999h . again pointing to a highly processed dust 
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Table 1. Characterizing parameters of HD 142527. ° | I I > 

Parameter Value 



Right Ascension 


IS" 56™4P:89 


Declination 


-42° 19' 23'.' 5 


Spectral Type 


F6IIIe 


Teff [K] 


6250 


Distance [pc] 


145±15 


Ay 


0.60±0.05 


Luminosity [Lq] 


20+2 


Radius [Rq] 


3.8±0.3 


Mass [Mo] 


2.2±0.3 


log Age [yr] 


6.7+0.4 


Group 


la 



Notes. We obtained new stellar parameters from our disk modeling ef- 
fort, see the text for details. 

environment. Recently. iHonda et alj (12009') showed evidence for 
crystalli ne water ice in the near-IR scattered light spectrum of 
the disk. IOhashil(l2008h showed that the cold dust is distributed in 
an arc-like structure, which adds to the complexity of the source. 

Clearly, the disk of HD 142527 is in an interesting evolution- 
ary phase. Its unusual properties might be caused by giant planet 
formation in the outer disk. In particular we are interested in 
understanding the very large far-lR excess, the geometry of the 
disk and the mineralogy of the dust. For this purpose we gath- 
ered a comprehensive data set on this object, consisting of SEST 
and ATCA millimeter photometry, a Spitzer-IRS spectrum, ISO- 
SWS and LWS spectra, VISIR N and Q-band imaging, VISIR 
N and Q-band spectroscopy, and optical, infrared and millime- 
ter photometry. In this paper, we focus on the SED and the disk 
geometry. In a future paper (Min et al., in preparation) we will 
study the dust mineralogy. This paper is structured as follows: 
In Sect. 12] we motivate the adopted stellar parameters. Section[3] 
describes the observations and the data reduction. In Sects. |4] 
through |6] we carry out an in-depth analysis of these data-sets. 
In Sect.|2]we construct a comprehensive model of the disk ge- 
ometry using a Monte Carlo radiative transfer code and address 
the discrepancies between model and observations. The discus- 
sion in Sect.[8]considers the emerging geometrical picture of the 
circumstellar matter and possible explanations for the peculiari- 
ties. Section|9] summarizes our main conclusions. 




30 20 10 -10 -20 -30 

RA offset (arosec; J20C0) 



Figure 1. ATCA 3.5 mm image centered on the peak of the emis- 
sion of HD 142527. The displayed contour levels are -16.2 (dot- 
ted line),16.2, 21.6, 27.0, 32.4, 37.8 mJy/beam. The resolution 
of the image is 16.2" by 2.9" at a PA of 24°. The beam size 
is shown in the bottom left. The source appears to be slightly 
resolved. 



Table 2. ATCA fluxes and positions. 



Name 


Peak flux [mjy] 




R.A. 




<5 




Uranus 


7.04 












IRSV 1540 


10.8 


15h 


44'" 39 !90 


-54° 


23' 


04'.'4 


HD 142527 


43.1±5.4 


15" 


56'" 41 IS? 


-42° 


19' 


23'.' 5 



Notes. The reported measurements of the calibrators and the target are 
obtained at 3.5 mm. For the phase-calibrator the flux was derived from 
imaging the six channels centered on the maser line. 



of ISiess et al. I (l2000h . which resulted in a stellar mass of M = 
2.2+0.3 Mq and an age of IQ^-^^'^^yr. Table[T] summarizes the 
stellar parameters. 



2. Stellar parameters 

HD 142527 w as cataloged as an F6III star b v iHoukI (Il978h 
[Henize ( 1976) noted its emission line nature and Waelke ns et al] 
(1996) clas sified it as a Herbig star With Hipparcos mea- 
surements, Ivan den Ancker et al I (119981) placed it at a dis- 
tance of 200^^° pc, b u t we consider the association by 
I Acke & van den Ancked (l2004t) to the star-formation region 
Sco OB2- 2 to be more reliable and therefore adopt the distance 
of 145 pc dde Zeeuw et a l] [T999l) . whi ch is within 2 s igma of the 
Hipparcos distance. By comparing a iKuruc j (Il99lh model for 
the photosphere associated with the given spectral type to the 
literature photometry, we obtained an optical extinction of Ay 
= 0.60+0.05 and a stellar luminosity of L = 15±2L0. However, 
our modeling of the dust (see Sect. 17. 2b demonstrates the pres- 
ence of a gray extinction component. Taking this component into 
account, we derive a stellar luminosity of L = 20+2 Lq. A stellar 
radius of - 3.8+0.3 Rq was derived from the effective tem- 
perature and the luminosity. The position in the Hertzsprung- 
Russell diagram was compared with the PMS evolution tracks 



3. Observations and data reduction 

3.1. SEST 

HD 142527 was observed with the 15 m Swedish/ESO sub- 
mm Telescope (SEST) in August 2003. The 37-chan nel SEST 
imaging bolometer array (SIMBA. lNvman et al.ll200ll) was used 
to observe the continuum dust emission at 1.2 mm. The ob- 
servations were performed in fast mapping mode, and maps 
were made of 400" by 500" in azimuth and elevation respec- 
tively. The atmospheric transparency was monitored by per- 
forming regular sky-dip measurements, the absolute flux cali- 
bration was established by observing Uranus. In total, six maps 
of HD 142527 were made, and the source was clearly detected 
in all of them. We determine the flux of HD 142527 at 1 .2 mm to 
be 1. 12+0.02 Jy. 

3.2. ATCA 

Observations at 3.476mm were obtained with the Australia 
Telescope Compact Array (ATCA) on the night of June 11, 
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Table 3. Log of the VISIR imaging. Table 4. Log of the VISIR spectroscopy. 



Target 


Filter 


A,- 


Chop 


Time 


Airmass 


Seeing 


Int.time 






[/J in] 


["] 


[h:in] 




["]" 


[s] 


HD 139127 


SiC 


11.85 


9 


05:07 


1.07 


0.83 


110 


HD 142527 


SiC 


11.85 


10 


05:18 


1.07 


0.83 


354 


HD 186791 


SiC 


11.85 


9 


10:07 


1.22 


0.77 


110 


HD 142527 


Q2 


18.72 


8 


02:21 


1.05 


0.82 


680 


HD 142527 


Q2 


18.72 


8 


02:37 


1.06 


0.75 


680 


HD 139127 


Q2 


18.72 


9 


02:57 


1.09 


0.72 


338 



Notes. The entries are ordered according to observation time. The 
columns from left to right give the target name, the imaging filter, its 
central wavelength, the chopper throw, the time of observation, the air- 
mass and optical seeing, and the effective integration time. 



2002 in the east-west 352 m configuration using a bandwidth of 
2.3 mm (see table|2]i. Only three antennae out of the six in the 
full array were available (numbers 2, 3, and 4). The resolution 
was 16.2" by 2.9". Observations on the target source were inter- 
leaved with observations of the bright, nearby SiO maser source, 
IRSV 1540, as a phase-reference calibrator, with a cycle of two 
minutes on the phase-reference followed by five minutes on the 
target. Every half hour T^y^ observations were made and every 
hour a new pointing solution derived using the phase-reference. 
The data reduction was completed using the MIRIAD data re- 
duction package (Sault et al.ll 1995b . The absolute flux scale and 
bandpass were calibrated using an observation of Uranus. The 
flux of Uranus was determined to be 7.0 Jy. From the variance 
in flux over time of the sources IRSV 1540 and 1057-797, the 
flux scale appears to be constant to within 10%. The measured 
central position matches the optical position of HD 142525 to a 
sub-arcsecond precision. Fitting a Gaussian to the source shown 
in Fig.[T]suggests that the source is marginally resolved, but the 
visibility data are too noisy to obtain a size estimate. The de- 
rived peak flux is 43.1 +5.4 mJy with a spatially integrated flux 
of 47.1±6.5mJv. This last number was previously quoted by 
lAcke et alj (|2004|P 1 as unpubHshed. 

3.3. VISIR imaging 

During the nights of April 28 and June 30, 2005 N and Q-band 
imaging of HD 142527 was performed with VISIR, the VLT 
Imager and Spectrometer for mid-IR (Lagage et al. 2004). The 
observations were obtained in visitor mode under the VISIR 
GTO program on circumstellar disks. The chosen filters were 
the SiC and Q2, which have central wavelengths at 11.85 and 
18.72/im and half-band widths of 2.34 and 0.88/im respec- 
tively. Standard chopping and nodding was applied (at 0.25 and 
0.033 Hz) to suppress the dominating atmospheric background 
emission. All the imaging data were obtained using a pixel 
scale of 0.075". The observing conditions during both nights 
were good: the airmass was below 1.1 and the optical seeing 
was around 0.8". For the imaging in the SiC filter the stan- 
dard stars HD 139127 (12.75 Jy) and HD 186791 (57 Jy) were 
bracketing the science observation and were used to estimate 
the photometry and the associated errors. The achieved sensi- 
tivity was estimated to be lOmJy/lOcrlh averaged over the two 
reference stars. For the imaging in the Q2 filter the reference 
star, HD 139127 (4.95 Jy), was used for PSF and photometric 
purposes. It was selected from the database of VISIR standard 
stars based on the flux level and distance to the science tar- 
get on the sky. The achieved sensitivity was estimated to be 

' In that paper the wavelength of the ATCA data-point was incor- 
rectly stated as 2.9 mm. 



Target 


R 




Orient. 


Tiine 


Airmass 


Seeing 


Int.time 






[/im] 


n 


[h:m] 




["] 


[s] 


HD 142527 


LR 


8.8 





05:58 


1.05 


0.76 


644 


HD 142527 


LR 


9.8 





06:13 


1.05 


0.76 


644 


HD 139127 


LR 


8.8 





06:38 


1.07 


0.67 


91 


HD 139127 


LR 


9.8 





06:42 


1.07 


0.60 


91 


HD 139127 


LR 


11.4 





06:45 


1.07 


0.61 


91 


HD 142527 


LR 


11.4 





07:01 


1.07 


0.67 


644 


HD 142527 


LR 


12.4 





07:17 


1.09 


0.67 


517 


HD 139127 


LR 


12.4 





07:40 


1.14 


0.73 


97 


HD 142527 


MR 


18.18 


113 


03:10 


1.11 


1.19 


1846 


HD 142527 


MR 


17.82 


113 


03:52 


1.18 


1.40 


1846 


HD 148478 


MR 


18.18 





05:17 


1.25 


0.86 


159 


HD 148478 


MR 


17.82 





05:23 


1.27 


0.83 


159 



Notes. The entries are ordered according to observation time. The 
columns from left to right give the target name, the spectral resolution 
(low/medium), the central wavelength, the orientation of the slit (rotat- 
ing from north to west), the time of observation, the airmass and optical 
seeing, and the efl^'ective integration time. 



70mJy/10crlh in the Q2 filter Table[3] summarizes the observa- 
tional details. For the reduction of the data we used a dedicated 
pipeline, which corrects for various in strumental signatures (see 
iPantin et al.ll2008ll2009l:IPantinll2010h . 

3.4. VISIR spectroscopy 

Long-slit spectra were taken with VISIR in the N-band and the 
Q-band. For both observations standard parallel chopping and 
nodding with a chopper throw of 8" was applied to correct for 
the atmospheric background. In table|4]we display the observa- 
tional log. 

The N-band spectrum was observed in the low-resolution 
mode (R~300) on the night of April 27, 2005. The 0.4" sHt 
was employed in a standard north-south orientation. The entire 
lOyum feature was covered in four separate settings. The observ- 
ing conditions were favorable; the average airmass was ~1.1 and 
the optical seeing was 0.7". The standard star HD 139127 was 
observed immediately before or after the science measurements 
for the correction of the atmospheric absorption. The N-band 
spectrum of HD 142527 was observed as part of a VISIR GTO 
program on circumstellar disks, in which t he 10 //m feature of a 
sample of 17 Herbig Ae stars was studied (i Verhoeff^et al]|20IOl) . 

The Q-band spectrum was observed in the medium- 
resolution mode (R~3000) on the night of June 30, 2007. The 
0.4" slit was employed in an inclined orientation with PA = 67°. 
The 18 jjm forsterite feature was mapped in two runs centered on 
17.82yum and 18.18jum. The observing conditions were fair; the 
average airmass was ~1.1 and the optical seeing was 1.4". The 
standard star HD 148478 was observed directly after the science 
measurement for the correction of the atmospheric absorption. 

The first steps of the data reduction are identical to the VISIR 
image reduction (see Sect. l3.3b . The resulting reduced spectral 
images contain the two-dimensional spectra in counts/s. These 
images were corrected for slit curvature and optical distortions. 
The spec tra wer e extracted using the optimal extraction tech- 
nique of iHornd (fl986). The wavelength calibration was per- 
formed by correlating an ALTRAN/HITR AN model wi th the 
simultaneously measured sky spectra. See iPantinI (1201 Ol) for a 
more detailed description of the reduction process. 

The N-band and Q-band spectrum were corrected for atmo- 
spheric extinction using one standard star and a simple correc- 
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Table 5. Photometric fluxes. 



2.5 3.1 



Band 




F[Jy] 


Reference 


Johnson U 


0.36 


0.34+0.01 


M98 


Johnson B 


0.44 


1.12±0.01 


M98 


Johnson V 


0.55 


1.84±0.02 


M98 


Near-IR J 


1.23 


3.8±0.1 


M98 


Near-IR H 


1.65 


4.6±0.1 


M98 


Near-IR K 


2.22 


5.6±0.1 


M98 


Near-IR L 


3.77 


7.4±0.3 


M98 


Near-IR M 


4.78 


6.7±0.3 


M98 


SiC 


11.85 


8.8±0.6 


VIO 


IRAS 12 


11.80 


10+4 


B88 


Q2 


18.72 


14.6±1.5 


VIO 


IRAS 25 


24.38 


21±5 


B88 


IRAS 60 


58.61 


105±12 


B88 


IRAS 100 


100.9 


84.7+16 


B88 


CSO0.8 


800 


4.2+0.5 


W95 


SEST 1.2 


1200 


1.12+0.02 


VIO 


CSO 1.3 


1300 


1.19+0.03 


W95 


ATCA 3.5 


3476 


0.047+0.006 


VIO 



S [Jy/arcsec'] 



References. B88 = iBeichman etal.llI988l: M98 = iMalfaitet alJl 19981 : 
VIO = this paper; W95 = IWalker & ButneillI995l . 



tion for the airmass difl'erence. The formahsm is defined as 

S A,c 



(1) 



where S and S are the observed counts of the cahbration 
and the science measurements, an d are the airmas s of the 
measurements, and / j ^ and , are the Coh en et al.l(ll999h model 
of the calibrator and the resulting spectrum. 

3.5. Additional data 

A lo w-resolution Spitzer- IRS spectrum of HD 142527 was taken 
from 'Juhasz et al.l (1201 Oh . The ISO spectrum was taken from 
iDominik et all (l2003h . N-band interferometric visibilities ob- 
served with MIDI were taken from Ivan Boekel et al ] (I2004ah . 
Photometric data points were taken from the l iterature and 
are lis ted in table |5] The optical photometry of Malf ait et al.l 
(119981) were converted to the standard Johnson UBV s ystem 
using the transformations of iRufener & Nicolei J 19881). The 



Subaru/Comics 24.5 /vm image was taken from Fuiiwa ra et al 



(2006) a nd the TIMMI2 N-band spectrum from , van Boekel et al 
(I2004bl) . 



4. Analysis of the millimeter data 

The two measurements at 1.2 and 3.5 mm together with two 
more measurements at 0.8 and 1.3 mm from the literature (see 
table |5]l give us an opacity index p {k oc A^^). Assuming the 
four data-points are in the Rayleigh-Jeans limit of the SED 
(J3 = a - 2), we obtain /3 = 1.0 + 0.1. Co mpared with other 
HAeBe stars this value is higher than average. ' Natta et aP (l2007h 
find that computed over the wavelength interval 1 .3 -7 mm, the 
opac ity index y6 is < 1 for about 60% of the objects. lAcke et all 
(1200 4) determined the spectral index in observations from 0.35 
to 2.7 mm and found an average of < > = 0.6 for 13 group I 
sources. At mm wavelengths the SED is thus steeper than that 
of most HAeBe stars. Consi dering that /3 - 1.7 for the ISM 
dWeingartner & Draind I20()l1) and assuming an optically thin 
disk at millimeter wavelengths, we can conclude that the dust 




Figure 2. Observed 18.72 fim VISIR image. The color bar shows 
the surface brightness in Jy/arcsec^. Overplotted are the posi- 
tions of the VISIR and TIMMI2 spectroscopic slits. The VISIR 
N-band slit position is plotted with the green straight line, the 
VISIR Q-band slit position is plotted with the blue dotted line 
and the TIMMI2 N-band spectrum is plotted with the yellow 
dashed line. The Spitzer slit encompasses the entire image. 



around HD 142527 has evolved, but that the average grain-size 
is relatively small compared with most HAeBe stars. 

From the 1.2 mm observation, we ca n derive a dust mass 
using the formalism of iHildebrandl (Il983h for an optically thin 
cloud 

FyCp- Aap 
B(v,T) 3Qy 

with a the grain radius and p the grain density, Qy the emissiv- 
ity, and d the distance to the star Considering the emissivity per 
unit dust mass Ky - 3Qy/4ap, taking Ky (250 fim) = lOcm^g"' 
from Hildebrand (1983) and extrapolating with p - 1.0 we ob- 
tain *■,,(!. 2mm) = 2cm^g"\ which is commonly used in the 
literature and should be correct within a factor of 4 (see e.g. 
iBeckwith et al.l[T990h . We assume an average dust temperature 
of Tdust - 40 K, which is not a large source of uncertainty, be- 
cause in the Rayleigh-Jeans limit the derived mass scales in- 
versely proportional with the temperature. Using the distance of 
d - 145 pc, we calculate a total dust mass in small grains of Mdust 
- 1 ■ 10"^ Mq. This value is consistent with the model we present 
in Sect.E] 



5. VISIR imaging analysis 

The photometric analysis yields integrated fluxes for HD 142527 
of 8.8+0.6Jy and 14.6+1.5 Jy in the SiC (11.85 yum, spectral 
resolution R = 5) and Q2 (18.72//m, R - 20) filters respec- 
tively. Because the emission of the star is negligible at these 
wavelengths, these fluxes can be attributed to the circumstel- 
lar disk. We searched for extended emission components at both 
wavelengths with a central component subtraction technique at 
a sub-pixel (1/10) precision level. The PSF in both filters was 



4 



A.P. Verhoeff et al.: The complex circumstellar environment of HD 142527. 




Figure 3. Normalized mean radial surface brightness profile of 
the VISIR 18.72//m image (black) and the PSF (red). We also 
plotted the result of our modeling effort (green; see Sect.|7]l. 




Figure 4. Central component-subtracted VISIR 18.72jum im- 
age of HD 142527. The color bar shows the surface brightness 
with a cut-off at 3.1 Jy/arcsec^. The overplotted contours from 
the 24.5 //m Subaru image are at 0.5, 1.0,2.0,3.0,4.0,5.0, and 
5.3 Jy/arcsec^. 



derived from the observation of standard stars just before or af- 
ter HD 142527. In the SiC filter, no visible extension is detected. 
This constrains the angular extent of the central component to 
< 0.35" in FWHM (50 AU at 145 pc). In the Q-band, extended 
emission is clearly detected. A visual inspection of the observed 
image in Fig.|2] demonstrates that the emission is divided be- 
tween a compact central component and an extended region. We 
find that the emission in the central component is also marginally 
spatially resolved (see Fig. [3]). A quadratic subtraction of the PSF 
leads to a FWHM of the spatial emission profile of 0.21"±0.04, 
which corresponds to 30 AU at a distance of 145 pc. 




Distance [' 



Figures. Mean radial surface brightness profiles of the east 
(PA = 0-180°) and west (PA = 180-360°) sides of the VISIR 
18.72/im image (blue) and the COMICS 24.5 /urn image (red). 
The top panel shows the profile of the observed images and the 
bottom panel shows the profile after central component subtrac- 
tion. 



The extended component that is visible after central com- 
ponent subtraction (figure |4]l has an elliptical shape. An ellipse 
fitting of the outer boundary isophotes indicates a position 
angle of 1 10°+20° and an aspect ratio of 1.3 (i = 40° ±20° from 
face-on). The integrated flux in the residuals after central com- 
ponent subtraction is 2.9±0.5Jy and the integrated flux of the 
subtracted central source is 1 1.1 + 1.1 Jy. However, the method 
used here of separating the photometry of the central component 
and the extended region is not very accurate, because their flux 
was mixed in the observation process by the PSF. In Sect. 15. 21 
we present results of a more precise deconvolution method. 



5.1. Comparison with previous Q-band imaging 

We compared our VISIR 18.72yu m image w ith the 

SUBARU/COMICS 24.5 fim image of iFuiiwara et all (l2006l) . 
applying the same central component subtraction technique 
as described above. In Fig.|4] we overplotted the central 
component-subtracted COMICS image with white contours. It 
shows that the eastern emission peak shifts to the south when 
going from 18.72 to 24.5 //m. This cannot be explained with a 
disk inclination or with flaring of the disk and is therefore an 
indication of azimuthal dust opacity variation, which must be 
linked to an azimuthal density variations. We wiU return to this 
point in Sect.lSH 

We also analyzed the mean surface brightness as a function 
of distance from the star In Fig.|5]we plot these spatial profiles 
for our VISIR image and for the COMICS image. The position 



5 



A.P. VerhoeiF et al.: The complex circumstellar environment of HD 142527. 



0.0 1.3 2.6 3.9 5.2 6.5 



S [Jy/orcsec^ 




200 
PA [degrees] 

Figures. Mean brightness profile of the outer disk (0.5" < r < 
1.1") as a function of the position angle (PA; east from north). 
The result of the deconvolved 18.72/im image is plotted in blue 
and the result of our model image is plotted in green. 



Figure 6. Deconvolved VISIR Q-band (18.72 /im) image of 
HD 142527. To zoom in on the lower surface brightness we 
made a cut-off at 6.5 Jy/arcsec^. 




1 .0 
Distance [' 



Figure?. Mean radial brightness profile as a function of dis- 
tance to the star as retrieved from the deconvolved 18.72 yum 
image (black). The result from the model image convolved with 
a Gaussian with a FWHM of 0. 1" is overplotted in green. 



of the star was assumed to coincide with the locus of the emis- 
sion optimum, which was derived to an accuracy of a tenth of a 
VISIR pixel with a two dime nsional Gaussian fit. F ollowing the 
same analysis approach as in iFujiwara et al.l (l2006h . we treated 
the eastern and western components separately: the eastern pro- 
file was averaged over the PA - 0-180° and the western pro- 
file was averaged over the PA - 180-360°. The observed pro- 
files in the top panel show that the flux ratio of /^i8.72/jm/^^24.5/im 
substantially decreases from the inner to the outer region. The 
bottom panel shows that the COMICS image is more extended. 
These features are consistent with a radially decreasing dust tem- 
perature. 



5.2. Deconvolution 

A more precise method to separate the central and outer emission 
components is to deconvolve the image. We dec onvolved the Q2 
image with the MEM-multi-resolution method (IPantin & Starckl 
1996) and obtained the image shown in Fig.|6l The spatial res- 
olution in the deconvolved image was limited to 0. 1 ", which is 
on the order of the Shannon limit imposed by the spatial sam- 
pling. At that resolution the image is composed of two compo- 
nents. The central component is marginally resolved (FWHM of 
0.15") and has an integrated flux of 10.84 Jy. The second com- 
ponent, located farther away at a distance of ~0.8", has a narrow 
ring-like shape. The ridge drawn by the maximum of emission 
in the ring can be fitted by an ellipse with an aspect ratio of 
1.2 (/ = 20° ±10° from face-on) and a position angle of PA = 
120°+20°. These numbers are slightly different from those ob- 
tained from the central component-subtracted image. We will re- 
turn to this in Sect. l7.5l The emission ring has a mean distance 
to the star (measured on the semi-major axis) of 0.8" (130 AU) 
and a FWHM of 0.2" (30 AU, see Fig.©. Figure|7] displays the 
mean radial profile of the deconvolved image, obtained by aver- 
aging the emission over 360°. The surface brightness as a func- 
tion of azimuth averaged over the radii from 0.5" to 1.1" has 
minima of 0.6 Jy/arcsec^ at PA = 0°, 180°, and 230° and a max- 
ima of 2.4 Jy/arcsec^ at PA - ~6Q° and 1.7 Jy/arcsec^ at PA = 
~300°(see Fig.^. 



6. VISIR spectroscopy analysis 

6.1. N-band 

To achieve a continuous VISIR low-resolution N-band spectrum, 
we spliced the four wavelengths settings together. We scaled 
them with the ratio of the median values of the overlapping re- 
gions. The flux of the spectrum was then scaled to the flux of 
the Spitzer spectrum at 10.6//m with an average factor of 1.24. 
We present the resulting spectrum in Fig.|9] We also plotted the 
Spitzer and TIMMI2 spectra for comparison. The excess at the 
blue end of the Spitzer spectrum is suggestive of extended low 
surface brightness PAH emission. Although there is no direct 
observational evidence in our data, this extended emission could 
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Figure 9. VISIR N-band spectrum of HD 142527 (black). 
Overplotted are the Spitzer spectrum (thick green) and the 
TIMMI2 specti-um (thin blue). 



Figure 11. Spatial emission profile of the Q-band spectrum at 
17.82yum (AA - O.lSyum). Overplotted in green is the corre- 
sponding calibration measurement. 



10 



17.8 



18.0 



Figure 10. VISIR Q-band specti-um of HD 142527. Overplotted 
are the Spitzer spectrum (thick green) and the wavelengths of 
poor atmospheric transmission (yellow bars). 



be linked to the components visible at 18.72 and 24.5 //m that 
largely fall outside the VISIR and TIMMI2 slits (see Fig.lHl. The 
small differences at 9.2, 10.0 and at 12.6/im are not well under- 
stood. 



6.2. Q-band 

Figure|2]shows the chosen slit orientation of the VISIR medium- 
resolution Q-band spectrum. The spectrum contains a contribu- 
tion from the central component and a contribution from the 
eastern component. In Fig. [10] we plot the integrated spectrum. 
The spectral setting centered on 18.18/im has been scaled by a 
factor 1.15 to continuously connect to the spectral setting cen- 
tered on 17.82jum. The wavelength regions where the atmo- 
spheric transmission as measured with the standard star is less 
than 30% of the maximum are masked out with yellow bands. 
We overplotted in thick green the Spitzer spectrum, which shows 
a significant difference in slope compared with the VISIR spec- 
trum. The bluer slope of the VISIR spectrum is in accordance 
with expectations, because the VISIR slit has selected the inner. 



hotter parts of the disk, while the Spitzer spectrum covers the 
entire disk, including all colder outer parts. 

The chosen slit orientation allows us to spatially separate the 
spectra of the central and eastern components. This is illustrated 
in Fig.[TT] where we plotted the spatial profile of the spectral 
section centered on 17.82/vm after a collapse of the wavelength 
dimension. In green we overplotted the calibrator. The eastern 
(left) and western (right) components are both significantly de- 
tected. Even the central component seems to be marginally re- 
solved, although the statistics are poor. We also extracted from 
the 2D spectral image two five pixel wide traces centered on the 
central emission and the eastern peak respectively. Their ratio 
showed no reliable spectral structures. 



7. Modeling 

To better understand the circumstellar matter around HD 142527 
we created a detailed geometric model that is capable of repro- 
ducing the described observations. We begin by explaining the 
principles of the model, the basic assumptions, and the setup. 
Then we constrain the free parameters and construct the final 
best-fit model using radiative transfer calculations and the ob- 
servations. This final model is tuned to accurately fit the SED, 
the Spitzer and ISO IR spectra, and the general features of the 
IR images of the star. We conclude our modeling effort with a 
discussion of the discrepancies between the simulated and ob- 
served data. 



7. 1 . Model basics and setup 

Our model consists of an axisymmetric structure that surrounds 
a single star. This structure is defined by its density distribution, 
composition, inclination, position angle, and inner and outer ra- 
dius. We u se the Monte Carlo radiative transfer code MCMax 
dMin et al.l 1200 9) to obtain a corresponding temperature distri- 
bution and simulated observables. To constrain the model pa- 
rameters we start with a basic setup which is refined after the 
radiative transfer calculation and comparison of simulated data 
with observables. In the radiative transfer calculation we only 
take the flux from the central star into account. Possible heating 
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by viscous processes in the disk can be ignored because the mass 
accretion rate of the disk is negligible. 

For the basic setup of the surface density of the disk we use 
a simple power law 



30 



nr) = So • - 



(3) ^ 



where r denotes the distance to the central sta r, the exponent p 
is set to the commonly used value -1 (see e.g. lDullemond et aP 
l2006l) . and the scaling factor, So, determines the mass of the disk. 
It is clear from the images (Figs.|4]and|6]l that the surface density 
of the disk is not continuous. We therefore split the disk into an 
inner and an outer disk separated by a gap. The structure of the 
two regions is treated separately. 

The vertical density distribution of the disk is solved by iter- 
ating the density and temperature structures until they are self- 
con sistent under the assump tion of hydrostatic equilibrium (see 
e.g. iDuUemond et"ani2007t) . Very small grains couple well to 
the gas and thus follow this hydrostatic equilibrium density dis- 
tribution. However, larger grains may partly decouple from the 
gas and settle towards the midplane of the disk through gravi- 
tation. To simulate this we added an extra par ameter that sca les 
the height of the dust in the disk (^P; see Ratz ka et al.l 120071 : or 
lAcke et ani2009l) . For values of ^' lower than unity, this can 
be interpreted as a settling of the grains to the midplane, hence 
we refer to this as the sedimentation parameter 

We fixed the composition and the grain size and shape distri- 
bution of the silicate comp onent of the dust of the inner disk to be 
equal to that obtained bv Ivan Boekel et aTl (l2005l) . In that paper 
the continuum flux was represented by emission from a single 
blackbody, which is sufficient for modeling the 10 /im feature. 
Here, we use a continuum opacity in the form of amorphous car- 
bon grains with a variable abundance to improve the fitting accu- 
racy in the optical and millimeter regimes. For the shape of the 
carbon grains, we us ed the distribution of hollow spheres (DHS; 
see [Min et al2 2005h. and for the refractive index we used the 
data by Prei bisch et alJ (11993 ^. The continuum component could 
also be reproduced by small grains of metallic iron and/or iron 
sulfide, as well as large grains of various dust species. However, 
this will have no great influence on the modeled geometry. For 
the outer disk we used the same basic compositional setup ex- 
cept for the addition of small water-ice grains. This is motivated 
by the clear water-ice feature around 45 fim as seen in the ISO 
spectrum. 

We assume that all types of dust species in the disk have the 
same temperature as the gas. This can be taken to mean that the 
different species are in thermal contact. To simulate the pres- 
ence of composite grain-like ice-mantles on silicate and carbon 
grains we simply added the opacities of the individual species. 
It would of course be physically more realistic to take the inter- 
actions of the diff erent componen ts in the composite grain into 
account (see e.g. Min et al]|2008l) . which would likely lead to a 
slightly different best-fit composition and grain size distribution. 
However, for understanding the geometrical characteristics and 
the general composition of the disk, the simple approach suf- 
fices. 

The prime puzzle of HD 142527 is its extremely large flux in 
the far-IR. The wavelength-integrated flux in the IR regime is al- 
most equal to the integrated flux from the star (Fir - 0.92 ■ Ft). 
Any successful disk model will first and foremost need to ex- 
plain this remarkable observation. Secondly, the peculiar spatial 
distribution of the IR emission as seen at 18.72 and 24.5 fim will 
need to be reproduced and the mid-IR interferometric visibilities 
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Figure 12. Mean radial brightness profiles of the east (0° < 
PA < 180°) and west (180° < PA < 360°) sides of the VISIR 
18.72/7m image (blue) and the COMICS 24.5 /urn image (red) 
compared with the profiles of the model images (green). 



have to be matched as well. Below we discuss several options to 
correctly model the SED and show how we arrived at our final 
best-fit model using the IR images and mid-IR interferometric 
constraints. 



7.2. Structure of the outer disk 

The central ingredient to obtain a high Fir/F, ratio is to move 
the bulk of the disk mass a large radial distance from the star 
where it can form a very high vertica l wall and rep rocess a large 
fraction of the stellar flux (see e.g. lDerodl2b07h . This can be 
understood in the following way. The fraction of the solid an- 
gle covered by the disk is approximately H/r, where H is the 
height of the optically thick disk over the midplane and r is the 
distance from the star. H can be written as the pressure scale 
height hp times an approximate constant (IChiang & GoldreichI 
il997i) . The standard disk theory shows that hp - c^/Qk where 
Cs - -^JkT /(juifip) is the isothermal sound speed at temperature 

T, with jdirip being the mean molecular weig ht. Qk = ^/GMJ? 
is the local Kepler frequency. The temperature of a freely irradi- 
ated wall decreases as a function of distance with r^aii r"'^^, 
so Cj oc r"'/'* and H/r oc r'^'*. This means that the larger the 
distance from the star, the higher H/r and Fir/F* will be. 

We find that a gap extending to 130 AU will allow the disk 
wall to rise high enough to account for most of the IR excess. 
The inclination and position angle of the disk were constrained 
by fitting the IR images and their radial surface brightness pro- 
files (Figs. [3] H] |6] |7] and \T2[ . The arc-like structures are well 
reproduced by the radiation from the inner edge of an inclined 
outer disk (see Fig. [T3T l. The outer radius was constrained by 
fine-tuning the fit to the narrow bump around 40 fim in the SED 
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Table 6. The parameters derived with the best fit model. 



Model 18.72 



4V 



COMICS 24.5 



Model 24.5 



Ik 



Figure 13. Compai-ison of the VISIR 18.72//m (upper left), 
COMICS 24.5 yum (bottom left), and their modeled images (to 
the right). The model images were obtained by convolving the 
model output at the specified wavelengths with the observed 
PSFs. For a clear comparison the observed and modeled images 
were displayed with the same lower and upper flux level cut-offs. 
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Figure 14. Model fit (black) to the literature photometry 
(squares with error bars) and ISO spectrum (gray). The differ- 
ent flux components, indicated in the bottom left, represent the 
photons whose last contact was with that particular model com- 
ponent. 



(see Fig.[T4ll. The disk mass was constrained by fine-tuning the 
fit of the flux level of the SED fit in the millimeter regimes. We 
note that for the necessary scale height of the outer disk to be 
achieved in a natural way, it is important that the stellar radi- 
ation reaches it largely unhindered. As we will see, combined 
with other observations, this poses significant constraints on the 
structure of the dust close to the star. 



Parameter 


Inner disk 


Halo 


Outer disk 


Inner radius (^i„ [AU]) 


0.3 


0.3 


130 


Outer radius (^om [AU]) 


30 


30 


200 


Inclination (/) 


20° 


- 


20° 


Position Angle (PA) 


160° 


- 


160° 


Sedimentation par. (T) 


0.4 




0.9 


Power law exponent (p) 


-1 


-1 


-1 


Dust mass {M [Mq]) 


2.5-10-' 


1.3-10-'° 


1.0-10-3 


Silicate abundance [%] 


63 


63 


34 


Carbon abundance [%] 


37 


37 


20 


Ice abundance [%] 






45 


Grain size [^um] 


0.1,1.5* 


0.1,1.5* 


2.5 


Crystallinity [%] 


23 


23 






Notes. Abundances a re given as mass fractions. (*) Following 
Van Boekel et al.l ( l2005h . 



■ MIDI observation 
Disk+Halo 
Disk 
Halo 




Figure 15. MIDI visibilities compared with the different models 
of the inner disk region. The model that combines an optically 
thick disk with an optically thin spherical halo fits the visibilities 
best. Typical uncertainties on the MIDI visibilities are below 5%. 



7.3. Structure of the inner disk 

The inner radius is constrained by fine-tuning the SED fit in the 
near-IR regime. For the vertical distribution of the material close 
to the star we consider two scenarios. In the first the material is 
confined in an optically thick disk and in the second the material 
is confined in an optically thin dust component with very large 
scale height. 

An optically thick inner disk? When we include a standard 
inner disk in the radiative transfer, the disk is strongly optically 
thick and casts a shadow on the outer disk. Furthermore, we need 
to artificially increase the scale height of this disk by roughly a 
factor of 3. While this produces the coiTect near-IR flux, it also 
casts a wide shadow onto the outer disk. The resulting loss of 
irradiated outer disk surface, and consequently the decrease in 
temperature, decreases the scale height and emission from the 
outer disk. To preserve the match with the observed mid- and far- 
IR flux, we would need to artificially increase the scale height of 
the outer disk as well, and that by a factor of 2.5 over the self- 
consistent hydrostatic equilibrium value. A model like this has 
several shortcomings. First of all, there is no physically realistic 
explanation for the increased scale height of the outer disk over 
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hydrostatic equilibrium, because this would require very high 
gas temperatures in the mid-plane, for which we have no justi- 
fication. Second, the resulting model images are more axisym- 
metric than the observed IR ima ges. And third, we compared the 
mid-lR visibilities obtained by Ivan Boekel et al. I t004a ) with 
MIDI with those obtained from this model and found that the 
model visibilities are a factor 1.5 too high (see Fig.fTSll. This 
means that the emission is too much peaked toward the center of 

the disk. 

A halo-like component? As suggested by iMulders et al.l 
(1201 oh . a way to avoid the extended shadow on the outer disk is 
to put the material in the inner region into an optically thin dust 
compone nt, which extends to a very high altitude, for example 
in a halo (IVinkovic et al.l2003ll2006h . Although halos have been 
the topic of much controversy, there are also other stars that 
show some indication of optically thin material e xtending to a 
high altitude in the inner disk (e.g. HD 163209. iBenistv et al.l 
12010) . We will discuss the possible origin of these halo-like 
structures in Sect. 18. 21 For our radiative transfer simulations, 
we modeled this component as a spherical halo, but we do not 
mean to imply that it has to be fully spherical in reality. Using 
an inner dust shell we can obtain a very good fit to both the 
images and the SED. Because we kept the composition and 
grain size distribution of the halo the same as for the inner 
disk, the halo causes a measure of gray extinction at optical 
wavelengths (r = 0.37). This means a slightly higher stellar 
luminosity L = 20 Lq is required than would be the case when 
only taking the interstellar extinction into account. The resulting 
slight decrease in scale height of the outer disk by a factor of 
0.8 can be easily understood in terms of settling of the grains 
toward the midplane. However, in order to fit the spectral shape 
(i.e. the range of temperatures) of the emission from the inner 
region, the powerlaw for the surface density needs to be very 
flat (2(r) oc r""'*). This leads to mid-IR visibilities that are a 
factor 4 lower than the observed visibihties (see Fig.fTSll. 



7.4. Final model 

We saw that the optically thick inner disk or the halo-like com- 
ponent alone do not deliver a satisfactorily model. The solution 
is to combine a low-mass inner disk in hydrostatic equilibrium 
with an optically thin halo-like component of dust. This allows 
us to obtain a very good agreement for the SED, the IR images 
and the mid-IR visibilities. 

We summarize our final best-fit model in Fig.[T7]and table|6] 
We assume an axisymmetric disk that is slightly inclined. It has 
a settled inner disk, which extends from 0.3 to 30 AU, and a 
sedimentation parameter of "I* = 0.4. At 0.3 AU the disk has 
a vertical height of 7 • 10"^ AU above the midplane (as traced 
by the radial t-\ surface in the optical). The dust reaches a 
temperature of ~1600K, which we deem acceptable. There is 
a small mass in an optically thin spherical halo between 0.3 AU 
and 30 AU. However, the outer radius of this halo is not well 
constrained by the observations. The outer disk starts at 130 AU. 
Vertically it was scaled to *!* = 0.9 times the hydrostatic equi- 
librium height, which puts the vertical height of the disk wall at 
60 AU (t= 1 surface). The inner region contains silicate dust wit h 
the same composition as obtained by Ivan Boekel et al ] (I2005h . 
and we used 37% carbonaceous grains to represent all contin- 
uum opacity sources. In order to correctly reproduce the wave- 
length position of the slope of the SED in the far-IR, we set the 
size of the grains in the outer disk to 2.5 fim and used 20% car- 
bonaceous grains. We also added 45% of crystalline water-ice 
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Figure 16. Spitzer-IRS LR spectrum of HD 142527 (SNR ^ 
400). Indicated are the wavelengths of the PAH bands (at 6.2, 
7.9, 8.6, 11.2, 12.7, and 16.4 pm), and of the forsterite bands 
(at 11.3, 16.2, 19.5, 23.7, 27.6 and 33.6 //m). The photometric 
points (red) are from our VISIR imaging. In blue we overplotted 
the modeled spectrum. 



to the dust-mixture of the outer disk to reproduce the observed 
45 ;um ice feature. 

7.5. Comparison of the model with the observations 

The disk model described above is capable to generally repro- 
duce the various observations. We are confident that the small 
discrepancies between observational and modeling results reflect 
details in the geometry and/or mineralogy that will not strongly 
affect the emerging general picture of the system. Here we list 
the comparison and discrepancies. 

Photometry. The fit to the literature photometry and ISO 
spectrum is presented in Fig. [14] The flux level and general slope 
of the photometry, as well as the features of the ISO spectrum are 
reproduced. The model overestimates the long wavelength part 
of the ISO spectrum by about 20%. 

Spitzer. The fit to the Spitzer spectrum is presented in Fig. [16] 
The flux level and general shape are reproduced, but the IQpm 
feature is overestimated by about 10%. The forsterite features 
at longer wavelengths are underestimated because we did not 
include small crystalline grains in the outer disk. 

MIDI. In Fig. [15] the MIDI visibilities were compared with 
the different models. Our final model reproduces the spatial 
scales probed by these MIDI visibilities. The diff'erences are 
most likely caused by gradients in the crystallinity of the disk, 
which we did not take into account. 

VISIR and COMICS. The modeled VISIR 18.72yum and 
COMICS 24.5 fim images are compared with the observations in 
Fig. [13] The simulated images were obtained by taking the out- 
put of the model at the specified wavelengths and by convolving 
them with their associated PSF The overall spatial distribution 
of the emission matches very well. In Fig.[3]we zoom in on the 
radial profile of the inner disk as seen in the VISIR image. In 
Fig.[T2]we plot the radial profile of the mean surface brightness 
of the eastern and western sides of the images. There is a slight 
discrepancy in the flux ratio of the eastern lobe over the central 
component of the COMICS image. 
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Figure 17. Pictographic display of our model of the disk around 
HD 142527. 

VISIR deconvolved. FigureQ shows the radial profile of the 
mean surface brightness of the deconvolved VISIR image. To 
model the deconvolved image we convolved the model output 
at 18.72yL(m with a 2D Gaussian with a FWHM of 0.1". This 
way we mimicked the obtained resolution of the deconvolved 
image. The fit is quite excellent, although the emission in the 
gap is not reproduced. This likely indicates that some dust is 
present in the gap. Figure[8] shows the azimuthal profile of the 
outer disk that results from the deconvolved image. The PA and 
the brightness of the eastern lobe are quite well reproduced, but 
the brightness peak at the western side of the outer disk is not 
reproduced. Furthermore, the image analysis gives a larger PA 
of the disk itself then the model does. These are indications that 
the density distribution and geometry of the outer disk are more 
complex than assumed by our model (compare with Fig. [12]), we 
will return to this in Sect. l8.4l 

8. Discussion 

8.1. General picture 

Our modeling effort shows that the circumstellar matter around 
HD 142527 consists of three geometrically distinct components, 
as depicted in Fig.[T7] First, there is a small inner disk structure 
running from 0.3 to 30 AU, with a puffed up inner rim and a 
self-shadowed disk further out. Second, there is a small amount 
of dust in an optically thin component extending to extremely 
high altitude. The third component is a fairly massive outer disk, 
consisting of an extremely high wall at 130 AU, which is in hy- 
drostatic equilibrium, and a self-shadowed outer part running up 
to200AU. 

8.2. Inner disk and halo-like structure 

In our model the inner disk and the halo-like structure together 
are responsible for the near-IR flux, the slightly resolved inner 
component of the VISIR image (see Fig.O, most of the IQfim 
feature and the MIDI visibilities. The inner disk contains very 
little mass in small dust grains (Mdust = 2.5 ■ 10"^ Mq). One 
likely explanation is that all the material is already locked up 
in planetesimals or larger bodies. The degree of sedimentation 
of the dust compared with the gas (*? = 0.4) seems to point at 



an advanced state of evolution as well. Expressed in the radial 
T=l surface in the optical, the height of the material at 0.3 AU 
reaches up to 7 ■ 10"^^ AU. 

For the halo-like structure we estimated a dust mass of 
1.3 ■ 10 '" Mq, which is one order of magnitude less than the 
inner disk. However, the halo-dust particles dominate the near- 
IR flux. The origin of the halo-like structure is not clear. It is 
generally expected that a dust halo decays naturally from stellar 
radiation pressure and gravity, therefore a mechanism is needed 
to replenish it. A possible replenishment scenario can be con- 
structed with a dynamically highly excited debris disk. An at- 
tempt to model the effects of such a scenario is underway (Krijt 
& Dominik 2010, in prep.). An obser vational hint for this t ype of 
replenishment was already found by iMalfait et al.l ( Il998h . who 
observed variable obscuration. There is also the possibility that 
some dust is entrained in the upp er disk layers where gas and 
dust temperatures strongly differ jWoitke etal.l 120091) . and the 
vertical distribution of the gas is much more extended than in 
models like our s, where gas and dust are coupled everywhere 
dThi et alJl201 Cft. A third possibility may be a disk-wind carry- 
ing small dust grains. 

A side effect of the presence of the halo-like structure is the 
gray extinction at optical wavelengths and a resulting higher es- 
timate of the stellar luminosity as mentioned in Sect.|2] The val- 
ues of the stellar radius, mass, and age given in table[T]were also 
derived from this new luminosity. 

8.3. Disk gap 

Our analysis of the VISIR 18.72 fim image and the modeling ef- 
fort (see Fig.|7ll demonstrate the presence of an extensive disk 
gap ranging from 30 to 130 AU, which is relatively void of small 
dust grains. One might speculate about the presence of one or 
more large Jupiter-mass planets in the gap that are responsi- 
ble for the clearing of the material. In the light of rece n t find- 
ings this seems a ver y likely possibility. Marois et al. (2008|) 
and lKalas et alJ (l2008l) found Jupiter-like planet s in wide orbits 
around young A-type stars with direct imaging. iLagrange et al.l 
(|2009^ consider the presence of a Jupiter-like planet in the de- 
bris disk of /3 Pictoris very likely. Furthermore, models of the 
solar nebula generally rely on an early formation of Jupiter, 
which then promot es the formation of the terrestrial planets 
(lPollacketalJll996h . This argument favors the presence of plan- 
etesimals or larger bodies in the inner disk (see Sect. l8.2l ). 

To check whether photoevapora tion could b e responsi ble fo r 
the clearing of the gap we consulted lGorti & Hol lenbach ( 20091) . 
who describe the following scenario for the dispersal of a disk 
around a 1 PMS star In the outer regions of the disk FUV 
photoevaporation dominates accretion and the disk is expected 
to rapidly shrink to a truncation radius, r, ~ 150 AU. Under favor- 
able conditions a gap is created at 3-30 AU, after which the inner 
disk rapidly disappears on a viscous timescale of < 10^ yr The 
outer disk is subsequently rapidly eroded (< 10^ yr) from the in- 
side out. The authors estimate disk lifetimes of ~10'' yr Because 
the gap of HD 142527 is located at larger radii and there is still 
material in the inner disk, it is unlikely that photoevaporation is 
the dominant cause of the disk geometry. 

8.4. Outer disk 

The most distinct and curious aspect of the appearance of 
HD 1425 27 is its large I R exc ess (Fir = 0.92 • Ft), which was 
noted by IDominik et all (l2003h . We managed to self-consistently 
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explain this flux level by invoking the presence of an extremely 
high wall of material at 130AU which is in hydrostatic equi- 
librium. We found that the height of this material reaches up to 
60 AU, as expressed in the radial t-1 surface in the optical. This 
means that the outer disk covers 42% of the sky as seen from 
the star. It is this covering fraction that allows the disk to repro- 
cess so much stellar light and emit IR radiation in two relatively 
small solid angles (2x29%). This wall is also the cause of the 
bright spot in the IR images. A frontal irradiation over almost 
the entire vertical surface is required in order for this wall to 
puff up the way it does. This is guaranteed by the small solid 
angle of the flattened inner disk (2%). 

Further out the disk is self-shadowed and runs up to 200 AU. 
From our model we estimated the total dust mass in small grains 
in the outer disk to be 1.0 • 10"^ M©, consistent with our simple 
analysis of the millimeter data (see Sec.|4]i. Assuming a gas-to- 
dust ratio of 100, the total disk mass is well be low the gravita- 
tional instability limits (see e.g. lGammiell200ll) . but it is at the 
high end of the range in disk masses ( 1 '-^ to 10 ' M©) found for 
other Herbig stars bv lAcke et aP (l2004 . 

The dust of the outer disk is a lso responsible f or the scat- 
tering seen in the Subaru images of lFukagawa et al.l (2006). The 
larger outer radius stated by these authors is likely caused in part 
by PSF smearing and in part by a very low-density dust com- 
ponent. This component does not contribute to our observations, 
so we did not take it into account. The cold dust of the outer 
disk is very processed, because it contains both crystalline H2O 
ice (see also iMalfait etal.lll999D and forsterite. lOhashil (l2008h 
find convincing evidence for strong gas depletion on the basis 
of '^CO (3-2) emission, which is also an indication of advanced 
disk evolution. 

The overall spatial brightness distribution was fairly well 
modeled, but from the overlay of the images at 18.72 and 
24.5 yum in Fig.|4] it is clear that the peak of the emission of 
the eastern lobe shifts to the south when going to longer wave- 
lengths. This asymmetry cannot be explained with an axisym- 
metric density distribution; thus this is an indication of azimuthal 
variation in the density. Also, Figs. l6]and[T3] show that there are 
some discrepancies between the observed images and the mod- 
eled images in the azimuthal direction. Previous modeling ef- 
forts of protoplanetary d isks have made the existence of spiral 
density waves likely (see lDurisen et al.ll2007l for a review). It is 
possible that we are witnessing the signature of these waves here 
in HD 142527. 

To add to the complexity of the geometry, the sub-millimeter 
mapping by Ohashi (2008) showed that the spatial distribution 
of the bulk of the material is in an arc-like structure enclosing 
the central star Their image shows concentrations to the north- 
east and northwest and a surprising absence of material to the 
south. The similarity to our deconvolved 18.72/zm is striking 
(see Fig.|6]l. From our azimuthal analysis we concluded that the 
emission maxima are actually at PA = 60° and PA = -60° 
(see Fig.|8]l. The resemblance that comes to mind is that of the 
Trojan asteroids of Jupiter, who reside in the two Lagrangian 
points of stability that lie 60° ahead and behind the planet's or- 
bit. This would imply the presence of a large planetary body 0.8" 
(120 AU) north of the star Such a body could also cause the in- 
ner edge of the outer disk to have an eccentric shape. This could 
explain the discrepancy found in the PA of the disk between the 
modeling and the image analysis (Sect. l5.2] i. A modeling of the 
gravitational dynamics is called for to determine whether this 
configuration is physically possible and to constrain the exact 
orbital parameters and mass of such a body. 



9. Conclusions 

We presented an observational and modeling study of the cir- 
cumstellar environment of HD 142527. The main conclusions of 
our work can be summarized as follows: 

- We obtained observations with SEST and ATCA, which pin 
the millimeter flux to the optical position of HD 142527. 
From the spectral slope in the millimeter we derived an opac- 
ity index offi- 1.0+0.1, indicative of a relatively small grain 
size. From the SEST flux at 1 .2 mm we derived a total dust 
mass in small grains of 1.0- lO^^* M©. 

- Mid-infrared imaging of the disk surrounding HD 142527 
using VISIR confirms the presence of a large gap extend- 
ing from ~30 to ~130 AU. We find that the disk is unre- 
solved in the 10 //m window with no sign of emission from 
the outer disk. The emission in the VlSlR 18.72 yum image is 
dominated by the marginally resolved inner disk and by the 
inner rim of the outer disk. A comparison between VISIR 
18.72 yum and SUBARU 24.5 yum imaging dFuiiwara et al.l 
1^06) shows good agreement; the spatial scale of the emis- 
sion increases with wavelength, which is consistent with an 
outward decreasing dust temperature. 

- We modeled the VISIR, SUBARU, MIDI, Spitzer, ISO and 
millimeter data with the MCMax radiative transfer code 
dMin et al]|2009l) . The large infrared excess can be under- 
stood with a self-consistent disk model with the following 
properties: (i) a geometrically flat inner disk running from 
0.3 AU up to 30 AU, (ii) an optically thin halo-like compo- 
nent of dust in the inner disk regions, (iii) an extensive disk 
gap, and (iv) a massive outer disk running from 130 AU up 
to 200 AU. The inner rim of the outer disk is directly irradi- 
ated by the central star, which causes a huge scale-height (h 
= 60 AU). The derived mass in small dust grains is consistent 
with our millimeter analysis (1.0 • lO"-' M©). 

- As a consequence of the gray extinction caused by the halo- 
like component we obtained new stellar parameters, includ- 
ing a stellar luminosity of L = 20+2 L© and a stellar age of 
log t = 6.7+0.4 yr 

- The deconvolved VISIR image shows azimuthal intensity 
variations that suggest a more complex structure of the outer 
disk than can be accounted for with our model. This suggests 
a non-axisymmetric distribution of material. 

- The presence of the halo-like component, the low mass in 
small grains and the high degree of sedimentation in the in- 
ner disk suggest the presence of planetesimals or larger bod- 
ies in the inner disk. The large disk gap, the highly processed 
nature of the grain population and the large mass of the outer 
disk suggest that planet formation is also on-going in the 
outer disk. HD 142527 is thus likely to evolve into a system 
similar to those recently found aroun d young A-type stars 
dMarois et al ] l2008l:lKalas et al.ll2008h . 
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